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— j  Physical  simulation  of  a  series  of  four,  six  cubic  meter  Liquefied  Natural  Gas 
(LNG)  spills  on  water  was  provided  by  the  Meteorological  Wind  Tunnel  facilities  at 
Colorado  State  University.  Field  data  were  collected  from  spills  performed  at  Naval 
Weapons  Center,  China  Lake,  California,  in  Fall  1978.  The  simulation  test  series  was 
to  provide  field  test  planning  information,  extend  the  value  of  a  limited  set  of 
field  measurements,  and  validate  the  concept  of  physical  modeling  of  LNG  plume  dis¬ 
persion  as  a  predictive  hazard  analysis  tool. 

Two  test  series  were  conducted  to  aid  in  placement  of  field  concentration 
measurement  Instrumentation.  The  first  utilized  1:170  scale  model  of  the  China  Lake 
site  in  the  Meteorological  Wind  Tunnel  and  the  second  utilized  1:85  scale  model  in 
the  Environmental  Wind  Tunnel.  Dispersion  data  collected  In  the  four  field  tests  at 
China  Lake  were  extended  In  two  test  series  conducted  In  the  Environmental  Wind 
Tunnel  over  a  1:85  scale  model.  The  following  results  were  obtained: 

(1)  Comparison  between  dispersion  data  for  similar  test  parameters  but  at  two 
different  model  scales,  1:85  and  1:170  produces  similar  concentration  variation. 

(2)  The  China  Lake  surface  roughness  and  topography  causes  the  LNG  vapor  plume 
to  disperse  more  rapidly  than  would  be  experienced  over  flat  terrain. 

(3)  For  the  four  simulated  LNG  field  tests  the  arrival  time,  peak  concentration, 
maximum  peak,  etc.,  have  been  tabulated;  peak  ground  level  concentration  contours 
have  been  determined;  and  normalized  peak  concentration  coefficients  vs.  downwind 
distance  have  been  plotted. 
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EXECUTIVE  SU1-MARY 

As  part  of  a  program  to  evaluate  the  hazards  associated  with 
Liquid  Natural  Gas  (LNG)  spills  a  series  of  six  cubic  meter  LNG  spill 
tests  were  performed  at  China  Lake  Naval  Weapons  Center.  Simultaneously, 
physical  simulation  of  these  spills  were  provided  by  the  meteorological 
wind  tunnel  facilities  at  Colorado  State  University.  The  intent  of 
the  simulation  test  series  was  to  provide  field  test  planning  informa¬ 
tion,  to  extend  the  value  of  a  limited  set  of  field  measurements,  and 
to  validate  the  concept  of  physical  modeling  of  LNG  plume  dispersion  as 
a  predictive  hazard  analysis  tool. 

Two  test  series  were  conducted  during  the  preliminary  investigations 
to  guide  the  placement  of  field  instrumentation.  The  first  series 
utilized  a  1:170  scale  model  of  the  China  Lake  site  in  the  Colorado 
State  University  Meteorological  Wind  Tunnel;  whereas  the  second 
series  incorporated  a  1:85  scale  model  in  the  CSU  Environmental  Wind 
Tunnel.  Subsequently  two  test  series  were  performed  over  a  1:85 
scale  model  in  the  Environmental  Wind  Tunnel  to  extend  dispersion  data 
for  four  field  tests;  China  Lake  LNG-18,  LNG-19,  LNG-20,  and  LNG-21. 

These  test  programs  provide  the  following  results: 

1)  A  comparison  between  dispersion  data  for  similar  test 
parameters  but  at  two  different  model  scales,  1:85  and  1:170 
produces  similar  concentration  variation. 

2)  The  China  Lake  surface  roughness  and  topography  causes  the 

LNG  vapor  plume  to  disperse  more  rapidly  than  would  be  experien¬ 
ced  over  flat  terrain. 


3)  The  arrival  tine,  peak  concentration,  naximun  peak, 
etc.  have  been  tabulated  together  for  the  four  simulated 
LNG  field  tests. 

4)  Peak  ground  level  concentration  contours  are  provided 
for  the  four  simulated  LNG  field  tests. 

5)  Normalized  peak  concentration  coefficients  with  respect 
to  downwind  distance  are  plotted  for  the  four  simulated 
LNG  field  tests. 
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1.0  INTRODUCTION 


Natural  gas  is  a  highly  desirable  form  of  energy  for  consumption 
in  the  United  States.  Its  conversion  to  heat  energy  for  home  and 
industrial  use  is  achieved  with  very  little  environmental  impact,  and 
a  sophisticated  distribution  network  already  services  a  major  part  of 
the  country.  Recent  efforts  to  expand  this  nation's  natural  gas 
supply  include  the  transport  of  natural  gas  in  a  liquid  state  from 
distant  gas  fields.  Unfortunately  storage  and  transport  of  liquid 
natural  gas  may  include  a  relatively  large  environmental  risk 
(Fay,  1973;  Burgess,  1972).  To  transport  and  store  liquified  natural 
gas  (LNG)  it  is  cooled  to  a  temperature  of  -162°  C.  At  this  temperature 
if  a  storage  tank  on  a  ship  or  land  were  to  rupture  and  the  contents 
spill  out  onto  the  earth's  surface,  rapid  boiling  of  the  LNG  would 
ensue  and  the  liberation  of  a  potentially  flammable  vapor  would  result. 

It  is  envisioned  that  if  the  flow  from  a  rupture  in  a  full  LNG  storage 
tank  could  not  be  constrained  28  million  cubic  meters  of  LNG  would  be 
released  in  80  minutes  (AGA,  1974).  Past  studies  (Neff,  1976;  AGA,  1974) 
have  demonstrated  that  the  cold  LNG  vapor  plume  will  remain  negatively 
buoyant  for  a  majority  of  its  lifetime;  thus,  it  represents  an  extreme 
ground  level  hazard,  this  hazard  will  extend  downwind  until  the  atmos- 

•  •  •  •  •  e  • 

phere  has  diluted  the  LNG  vapor  below  the  lower  flammability  limit  (a 
local  concentration  for  methane  below  5%  by  volume) . 

It  is  important  that  accurate  predictive  models  for  LNG  vapor  cloud 
physics  be  developed,  so  that  the  associated  hazards  of  transportation 
and  storage  may  be  evaluated.  Various  industrial  and  governmental 
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agencies  have  sponsored  a  combination  of  analytical,  empirical, 
and  physical  modeling  studies  to  analyze  problems  associated  with  the 
transportation  and  storage  of  LNG.  Since  these  models  require  assumptions 
to  permit  tractable  solution  procedure  one  must  perform  atmospheric 
scale  tests  to  verify  their  accuracy. 

The  objectives  of  the  present  study  were  to  help  the  Coast  Guard 
and  other  contractors  prepare  for  a  series  of  atmospheric  scale  LNG 
spill  tests  at  the  China  Lake  Naval  Weapons  Center  site  and,  subsequently 
to  model  the  field  program  at  reduced  scales  in  meteorological  wind  tunnels. 
This  test  validation  program  was  designed  to  examine  the  validity  of 
physical  modeling  of  an  LNG  vapor  plume  and  to  document  any  limitations 
of  this  technique. 

The  physical  modeling  study  was  conducted  in  two  phases.  The  first 
phase  was  a  pre-test  series  which  provided  general  information  that  was 
needed  in  both  the  laboratory  and  field  to  sqt  up  instrumentation  and 
investigate  modeling  techniques.  The  second  phase  was  a  post-test 
series  in  which  four  different  spill  tests  were  simulated  in  the  wind 
tunnel.  The  meteorological  and  source  conditions  for  these  four  tests 
are  summarized  in  Table  1. 

The  methods  employed  in  the  physical  modeling  of  atmospheric  and 
plume  motion  are  discussed  in  Chapter  2.  The  details  of  model  construc¬ 
tion  and  experimental  measurements  are  described  in  Chapter  3.  Chapter  4 
discusses  the  test  program  and  results  obtained. 
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2.0  MODELING  OF  PLUME  DISPERSION 

To  obtain  a  predictive  model  for  a  specific  plume  dispersion 
problem  one  must  quantify  the  pertinent  physical  variables  and 
parameters  into  a  logical  expression  that  determines  their  interela- 
tionships.  This  task  is  achieved  implicitly  for  processes  occurring 
in  the  atmospheric  boundary  layer  by  the  formulation  of  the  equations 
of  conservation  of  mass,  momentum,  and  energy.  These  equations  with 
site  and -source  conditions  and  associated  constituitive  relations  are 
highly  descriptive  of  the  actual  physical  interrelationship  of  the 
various  independent  (space  and  time)  and  dependent  (velocity,  tempera* 
ture,  pressure,  density,  etc.)  variables. 

These  generalized  conservation  statements  subject  to  the  typical 
boundary  conditions  of  atmospheric  flow  are  too  complex  to  be  solved 
by  present  analytical  or  numerical  techniques.  It  is  also  unlikely  that 
one  could  create  a  physical  model  for  which  exact  similarity  exists 
for  all  the  dependent  variables  over  all  the  scales  of  motion  present 
in  the  atmosphere.  Thus,  one  must  resort  to  various  degrees  of  approxi¬ 
mation  to  obtain  a  predictive  model.  At  present  purely  analytical  or 
numerical  solutions  of  plume  dispersion  are  unavailable  because  of  the 
classical  problem  of  turbulent  closure  (Hinze,  1975).  Such  techniques 
rely  heavily  upon  empirical  input  from  observed  or  physically  modeled 
data.  The  combined  empirical-analytical-numerical  solutions  have  been 
combined  into  several  different  predictive  approaches  by  Pasquill  (1974) 
and  others.  The  estimates  of  dispersion  by  these  approaches  are  often 
crude;  hence,  they  should  only  be  used  when  the  approach  and  site  terrain 
are  uniform  and  without  obstacles.  Boundary  layer  wind  tunnels  are 


capable  of  physically  modeling  plume  processes  in  the  atmosphere 
under  certain  restrictions.  These  restrictions  are  discussed  in  the 


next  few  sections. 


2.1  Physical  Modeling  of  the  Atmospheric  Boundary  Layer 

The  atmospheric  boundary  layer  is  that  portion  of  the  atmosphere 
extending  from  ground  level  to  approximately  100  meters  within 
which  the  major  exchanges  of  mass,  momentum,  and  heat  occur.  This 
region  of  the  atmosphere  is  described  mathematically  by  statements 
of  conservation  of  mass,  momentum,  and  energy  (Cermak,  1971). 

The  general  requirements  for  laboratory-atmospheric-flow  similarity 
may  be  obtained  by  fractional  analysis  of  these  governing  equations 
(Kline,  1965).  This  methodology  is  accomplished  by  scaling  the 
pertinent  dependent  and  independent  variables  and  then  casting 
the  equations  into  dimensionless  form  by  dividing  through  by  one  of 
the  coefficients  (the  inertial  terms  in  this  case) .  Performing  these 
operations  on  such  dimensional  equations  yields  dimensionless  parameters 
commonly  known  as: 


Reynolds  number  Re  =  UqLo/vo 


Inertial  Force 
Viscous  Force 


Bulk  Richardson 
number 


Rossby  number 


Prandtl  number 


“  ■  [<STVTo]  (Lo/Uo>  *0 


Gravitational  Force 
Inertial  Force 


Ro 


U  /L  n 

O  0  0 


Inertial  Force 
Coriolis  Force 


■  W'oS  ’ 

o 


Viscous  Diffusivity 
Thermal  Diffusivity 


2  _ 

Eckert  number  Ec  -  U  /C  (AT) 
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For  exact  similarity  between  different  flows  which  are  described 
by  the  same  set  of  equations,  each  of  these  dimensionless  parameters 
must  be  equal  for  both  flow  systems.  In  addition  to  this  requirement, 
there  must  be  similarity  between  the  surface-boundary  conditions. 

Surface-boundary  condition  similarity  requires  equivalence  of  the 
following  features: 

a.  Surface-roughness  distributions, 

b.  Topographic  relief,  and 

c.  Surface-temperature  distribution. 

If  all  the  foregoing  requirements  are  met  simultaneously,  all 
atmospheric  scales  of  motion  ranging  from  micro  to  mesoscale  could 
be  simulated  within  the  same  flow  field  for  a  given  set  of  boundary 
conditions  (Cermak,  197s).  However,  all  of  the  requirements  cannot  be 
satisfied  simultaneously  by  existing  laboratory  facilities;  thus,  a 
partial  or  approximate  simulation  must  be  used.  This  limitation  requires 
that  atmospheric  simulation  for  a  particular  wind-engineering  application 
must  be  designed  to  simulate  most  accurately  those  scales  of  motion 
which  are  of  greatest  significance  for  the  given  application. 

2.1.1  Partial  Simulation  of  the  Atmospheric  Boundary  Layer 

A  partial  simulation  is  practically  realizable  only  because  the 
kinematics  and  dynamics  of  flow  systems  above  a  certain  minimum 
Reynolds  number  are  independent  of  this  number's  magnitude  (Schlichting, 
1968;  Zoric,  1972).  The  magnitude  of  the  minimum  Reynolds  number 
will  depend  upon  the  geometry  of  the  flow  system  being  studied. 

Halitsky  (1969)  reported  that  for  concentration  measurements  on  a  cube 


placed  in  a  near  uniform  flow  field  the  Reynolds  number  required  for 
invariance  of  the  concentration  distribution  over  the  cube  surface  and 
downwind  must  exceed  11,000.  Because  of  this  invariance  exact  similarity 
of  Reynolds  parameter  is  neglected  when  physically  modeling  the 
atmosphere . 

When  the  flow  scale  being  modeled  is  small  enough  such  that  the 
turning  of  the  mean  wind  directions  with  heights  is  unimportant, 
similarity  of  the  Rossby  number  may  be  relaxed.  For  the  case  of  dispersion 
of  LNG  near  the  ground  level  the  Coriolis  effect  on  the  plume  motion 
would  be  extremely  small. 

T 

The  Eckert  number  for  air  is  equivalent  to  0.4  Kr  (^r-) 

o 

where  Ma  is  the  mach  number  (Hinze,  1975).  For  the  wind  velocities 

and  temperature  differences  which  occur  in  either  the  atmosphere  or 
the  laboratory  flow  the  Eckert  number  is  very  small;  thus,  the  effects 
of  energy  dissipation  with  respect  to  the  convection  of  energy  is 
negligible  for  both  model  and  prototype.  Eckert  number  equality  is 
relaxed. 

Prandtl  number  equality  is  easily  obtained  since  it  is  dependent 
on  the  molecular  properties  of  the  working  fluid  which  is  air  for  both 
model  and  prototype. 

Bulk  Richardson  number  equality  may  be  obtained  in  special  labora¬ 
tory  facilities  such  as  the  Meteorological  Wind  Tunnel  at  Colorado 
State  University  (Plate,  1963). 
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Quite  often  during  the  modeling  of  a  specific  flow  phenomena 
it  is  sufficient  to  model  only  a  portion  of  the  boundary  layer  or  a 
portion  of  the  spectral  energy  distribution.  This  relaxation  allows 
more  flexibility  in  the  choice  of  the  length  scale  that  is  to  be  used 
in  a  model  study.  When  this  technique  is  employed  it  is  common  to  scale 
the  flow  by  any  combination  of  the  following  length  scales,  6,  the 
portion  of  the  boundary  layer  to  be  simulated;  zq,  the  aerodynamic 
roughness;  A^,  the  integral  length  scale  of  the  velocity  fluctuations, 
or  Ap,  the  wave  length  at  which  the  peak  spectral  energy  is  observed. 

Unfortunately  many  of  the  scaling  parameters  and  characteristic 
profiles  are  difficult  to  obtain  in  the  atmosphere.  They  are  infre¬ 
quently  known  for  many  of  the  sites  to  which  a  model  study  is  to  be 
performed.  To  help  alleviate  this  problem  Counihan  (1975)  has  summarized 
measured  values  of  some  of  these  different  discretions  for  the 
atmospheric  flow  at  many  different  sites  and  flow  conditions. 

2.2  Physical  Modeling  of  Plume  Motion 

In  addition  to  modeling  the  turbulent  structure  of  the  atmosphere 
in  the  vicinity  of  a  test  site  it  is  necessary  to  properly  scale  the 
plume  source  conditions.  One  approach  would  be  to  follow  the  methodo¬ 
logy  used  in  Section  2.1.,  i.e.,  writing  the  conservation  statements 
for  the  combined  flow  system  followed  by  fractional  analysis  to  find 
the  governing  parameters.  An  alternative  approach,  the  one  which  will 
be  used  here,  is  that  of  similitude  (Kline,  1965).  The  method  of 
similitude  obtains  scaling  parameters  by  reasoning  that  the  mass  ratios, 
force  ratios,  energy  ratios,  and  property  ratios  should  be  equal  for 
both  model  and  prototype.  When  one  considers  the  dynamics  of  gaseous 


plume  behavior  the  following  nondimensional  parameters  of  importance 
are  identified  (Hoot,  1974;  Skinner,  1978;  Snyder,  1972;  Halitsky,  1969) 
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to  obtain  simultaneous  simulation  of  these  four  parameters  it  is 
necessary  to  maintain  equality  of  the  plume's  specific  gravity 

ps/pa‘ 


It  has  been  assumed  that  the  dominant  transfer  mechanism  is  that  of 
turbulent  entrainment.  Thus  the  transfer  processes  of  heat  conduction, 
convection,  and  radiation  are  negligible. 

2 

The  scaling  of  plume  Reynolds  number  is  also  a  significant  parameter. 
Its  effects  are  invariant  over  a  large  range  thus  making  it  possible 
to  scale  the  distribution  of  mean  and  turbulent  velocities  and  relax 
exact  parameter  equality. 
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2.2.1  Partial  Simulation  of  Plume  Motion 

Thd  restriction  to  an  exact  variation  of  the  density  ratio  for 

the  entire  life  of  a  plume  is  difficult  to  meet  for  plumes  which 

simultaneously  vary  in  molecular  weight  and  temperature.  To  emphasize 

this  point  more  clearly,  consider  the  mixing  of  two  volumes  of  gas,  one 

being  the  source  gas,  ¥  ,  the  other  being  ambient  air,  ¥  .  Conside- 

s  a 

ration  of  the  conservation  of  mass  and  energy  for  this  system  yields 
(Skinner,  1978) : 1 


P. 


(2-7) 


If  the  temperature  of  the  air,  T  ,  equals  the  tejqperature  of  the  source 

a 

gases,  T  ,  or  if  the  product,  C  M,  is  equal  for  both  source  gas  and 
s  p 

air  then  the  equation  reduces  to: 


—  ¥  +  ¥ 
Pa  s  a 

¥  +  ¥ 
s  a 


(2-8) 


Thus  for  two  prototype  cases:  1)  an  isothermal  plume  and  2)  a  thermal 
plume  which  is  composed  of  air,  it  does  not  matter  how  one  models 
the  density  ratio  as  long  as  the  initial  density  ratio  value  is  equal 
for  both  model  and  prototype. 


*The  pertinent  assumption  in  this  derivation  is  that  the  gases  are 
ideal  and  properties  are  constant. 
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For  a  plume  whose  temperature,  molecular  weight,  and  specific 
heat  are  all  different  from  that  of  the  ambient  air,  i.e.,  a  cold 
natural  gas  plume,  equality  in  the  variation  of  the  density  ratio  upon 
mixing  must  be  relaxed  slightly  if  one  is  to  model  utilizing  a  gas 
different  from  that  of  the  prototype.*  In  most  situations  this  deviation 
from  exact  similarity  is  very  small.  (See  discussion  Section  3.2 
and  Figure  2.) 

Scaling  of  the  effects  of  heat  transfer  by  conduction,  convection, 
radiation,  or  latent  heat  release  from  entrained  water  vapor  cannot  be 
reproduced  when  the  model  source  gas  and  environment  are  isothermal. 
Fortunately  in  a  large  majority  of  industrial  plumes  the  effects  of  heat 
transfer  by  conduction,  convection,  and  radiation  from  the  environ¬ 
ment  are  small  enough  that  the  plume  buoyancy  essentially  remains 
unchanged.  The  influence  of  latent  heat  release  by  moisture  upon  the 
plume's  buoyancy  is  a  function  of  the  quantity  of  water  vapor  present 
in  the  plume  and  the  humidity  of  the  ambient  atmosphere.  Such  phase 
change  effects  on  plume  buoyancy  can  be  very  pronounced  in  some  prototype 
situations.  Figure  1  displays  the  variation  of  specific  gravity  from 
a  spill  of  liquified  natural  gas  in  atmospheres  of  different  humidities. 

A  reasonably  complete  simulation  may  be  obtained  in  some  situations 
even  when  modified  density  ratio  p  /p  is  stipulated.  The  advantage 
of  such  a  procedure  is  demonstrated  most  clearly  by  the  statement  of 
equality  of  Froude  numbers. 


If  one  was  to  use  a  gas  whose  temperature  is  different  from  that  of 
the  ambient  air  then  consideration  of  similarity  in  the  scaling  of  the 
energy  ratios  must  be  considered. 
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Solving  this  equation  to  find  the  relationship  between  model  velocity 
and  prototype  velocity  yields: 


where  S.G.  is  the  specific  gravity,  CP  /p  )»  and  L.S.  is  the  length 
scale,  (Lp/L^) .  By  increasing  the  specific  gravity  of  the  model  gas 
compared  to  that  of  the  prototype  gas,  for  a  given  length  scale,  one 
increases  the  reference  velocity  used  in  the  model.  It  is  difficult 
to  generate  a  flow  which  is  similar  to  that  of  the  atmospheric  boundary 
layer  in  a  wind  tunnel  run  at  very  low  wind  speeds.  Thus  the  effect  of 
modifying  the  models  specific  gravity  extends  the  range  of  flow  situations 
which  can  be  modeled  accurately.  But  unfortunately  during  such  adjust¬ 
ment  of  the  model  gas's  specific  gravity  at  least  two  of  the  four 
similarity  parameters  listed  must  be  neglected.  The  options  as  to  which 
two  of  these  parameters  to  retain,  if  any,  depends  upon  the  physical 
situation  being  modeled.  Two  of  the  three  possible  options  are  lifted 
below. 


(1)  Froude  No.  Equality 

Momentum.  Ratio  Equality 
Mass  Ratio  Inequality  ^ 
Velocity  Ratio  Inequality 


*When  this  technique  is  employed  distortion  in  velocity  scales  or 
similarly  volume  flow  rates  requires  a  correction  in  source  strength. 
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(2)  Froude  No.  Equality 

Momentum  Ratio  Inequality 
Mass  Ratio  Inequality 
Velocity  Ratio  Equality 

Both  of  these  schemes  have  been  used  to  model  plume  dispersion  downwind 
of  an  electric  power  plant  complex  (Skinner,  1978)  and  (Meroney,  1974) 
respectively. 

The  modeling  of  the  plume  Reynolds  number  is  relaxed  in  all 
physical  model  studies.  This  parameter  is  thought  to  be  of  small 
importance  since  the  plume’s  character  will  be  dominated  by  background 
atmospheric  turbulence  soon  after  its  emission.  But,  if  one  was 
interested  in  plume  behavior  near  the  source,  then  steps  should  be  taken 
to  assure  that  the  model's  plume  is  fully  turbulent. 

2.3  Modeling  of  Plume  Dispersion  at  China  Lake 

In  the  sections  above  a  review  of  the  extent  to  which  wind  tunnels 
can  model  plume  dispersion  in  the  atmospheric  boundary  layer  has  been 
presented.  In  this  section  these  arguments  will  be  applied  to  the 
specific  case  of  an  LNG  spill  at  the  China  Lake  Naval  Weapons  Center. 

2.3.1  Physical  Modeling  of  the  China  Lake  Atmospheric  Surface  Layer 
Only  small  quantities  of  LNG  are  being  spilt  at  the  China  Lake 
test  site,  approximately  five  cubic  meters;  hence,  it  was  decided  to 
scale  the  topography  at  the  site  by  a  one-to-eightyfive  ratio  to  improve 
the  resolution  during  concentration  measurements.  At  this  scale  it  is 
not  possible  to  simulate  the  entire  depth  of  the  atmospheric  boundary 
layer.  This  lack  of  simulation  of  the  entire  boundary  layer  is  consi¬ 
dered  to  be  of  minor  consequence  since  the  source  is  at  the  ground  level 
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and  the  vertical  dispersion  in  the  near  field  is  minimal.  The 
aerodynamic  roughness  (zq),  power  exponent  (1/n),  and  the  integral 
length  scales  of  turbulence  were  all  scaled  to  prototype  equivalent 
values  at  a  two  meter  height.  The  specific  values  obtained  are 
listed  in  Table  6.  A  preliminary  validation  program  was  conducts  at 
a  model  scale  of  1:170  to  assure  that  this  partial  depth  simulation 
method  yielded  reasonably  consistent  results  over  two  different 
scale  ratios. 

The  equality  of  Richardson  numbers  was  not  specified.  During 
tests  simulated  herein  the  prevailing  prototype  conditions  were  either 
neutral  or  only  slightly  unstable;  thus,  buoyancy  generated  turbulence 
played  only  a  small  role  in  determining  ground  level  turbulent  structure. 

2.3.2  Physical  Modeling  of  the  China  Lake  LNG  Spill  Plume 

The  buoyancy  of  a  plume  resulting  from  an  LNG  spill  is  a  function 
of  both  the  mole  fraction  of  methane  and  temperature.  If  the  plume 
entrains  air  adiabatically,  then  the  plume  would  remain  negatively 
buoyant  for  its  entire  lifetime.  If  the  humidity  of  the  atmosphere  were 
high  then  the  state  of  buoyancy  of  the  plume  will  vary  from  negative 
to  weakly  positive.  These  conclusions  are  bom  out  in  Figure  1,  which 
illustrates  the  specific  gravity  of  a  mixture  of  methane  at  boiloff 
temperature  with  ambient  air  and  water  vapor. 

Since  the  adiabatic  plume  assumption  will  yield  the  most  conservative 
downwind  dispersion  estimates  this  situation  was  simulated.  Several 
investigators  have  confirmed  that  the  Froude  number  is  the  parameter 
which  governs  plume  spread  rate,  trajectory,  plume  size  and  entrainment 
during  initial  dense  plume  dilution  (Hoot  and  Meroney,  1974;  Bodurtha, 
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1961;  Van  Ulden,  1974;  Boyle  and  Kneebone,  1973).  The  modeling  of 
momentum  is  not  of  critical  importance  for  a  ground  source  released 
over  a  fairly  large  area.  The  equality  of  model  and  prototype  specific 
gravity  was  relaxed  so  that  pure  Argon  gas  could  be  used  for  the  model 
source  gas. 

Argon  provides  almost  eight  times  the  detection  sensitivity  for 
instantaneous  concentration  measurements  as  the  carbon  dioxide  used  in 
previous  studies  (Meroney,  1977).  The  variation  of  specific  gravity 
with  equivalent  observed  mole  fraction  of  methane  for  these  different 
gases  is  plotted  in  Figure  2.  The  variation  of  Proude  number  with 
equivalent  mole  fraction  of  methane  for  the  simulation  gas  used.  Argon, 
is  plotted  in  Figure  3.  Over  the  concentration  range  where  the  buoyancy 
forces  are  dominant  the  variation  of  the  Froude  number  is  properly 
simulated.  Undistorted  scaling  of  velocity  components  was  maintained, 
which  implies  the  undistorted  scaling  of  source  strength. 

The  actual  source  condition,  boiloff  rate  per  unit  area  over  the 
time  duration  of  the  spill,  for  a  spill  of  LNG  on  water  is  highly 
unpredictable.  As  there  was  no  data  on  the  variable  area  and  variable 
volume  nature  of  the  different  LNG  tests  conducted  at  China  Lake  the 
source  conditions  were  approximated  by  assuming  a  steady  boiloff  rate 
for  the  duration  of  the  spill  over  a  constant  area. 

Since  the  thermally  variable  prototype  gas  was  simulated  by  an 
isothermal  simulation  gas,  the  concentration  measurements  observed  in 
the  model  must  be  adjusted  to  equivalent  concentrations  that  would  be 
measured  in  the  field.  This  relationship  which  is  derived  in  Appendix  A 
is : 
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\  ’  (1  -  Xm'r 

a 


where 

X  =  volume  or  mole  fraction  measured  during  the  model  tests 
m 

Tg  =  source  temperature  of  LNG  during  field  conditions 
and  =  ambient  air  temperature  during  field  conditions 


16 


3.0  DATA  ACQUISITION  AND  ANALYSIS 

The  methods  used  to  make  laboratory  measurements  and  the  techniques 
used  to  convert  these  measured  quantities  to  meaningful  field  equivalent 
quantities  are  discussed  in  this  section.  Attention  has  been  drawn  to 
the  limitations  in  the  techniques  in  an  attempt  to  prevent  misinterpreta¬ 
tion  or  misunderstanding  of  the  results  presented  in  the  next  section. 
Some  of  the  methods  used  are  conventional  and  need  little  elaboration. 

3.1  Wind  Tunnel  Facilities 

The  first  test  series  of  the  experiments  were  performed  in  the 
Meteorological  Wind  Tunnel  (MWT)  shown  in  Figure  4.  This  wind  tunnel, 
especially  designed  to  study  atmospheric  flow  phenomena,  incorporates 
special  features  such  as  an  adjustable  ceiling,  temperature  controlled 
boundary  walls,  and  a  long  test  section  to  permit  reproduction  of 
micrometeorological  behavior.  Mean  wind  speeds  of  0.3  to  40  m/s 
can  be  obtained  in  the  MWT.  Boundary- layer  thickness  up  to  1.2  m  can 
be  developed  "naturally"  over  the  downstream  6  m  of  the  MWT  test  section. 
Thermal  stratification  in  the  MWT  is  provided  by  the  heating  and 
cooling  systems  in  the  section  passage  and  test  section  floor.  The 
flexible  test  section  roof  on  the  MWT  is  adjustable  in  height  to  permit 
the  longitudinal  pressure  gradient  to  be  set  at  zero.  The  MWT  facility 
is  described  in  detail  by  Plate  and  Cermak  (1963). 

A  set  of  vortex  generators  was  installed  0.6  m  downwind  of  the 
entrance  to  give  the  simulated  boundary  layer  an  initial  impulse  of 
growth.  These  vortex  generators  were  then  followed  by  9  m  of  smooth 
floor  and  a  3  m  approach  ramp  to  the  1:170  scaled  topography  of  the 


China  Lake  site. 
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The  Environmental  Wind  Tunnel  (EWT)  shown  in  Figure  5  was  used  for 
the  remaining  three  test  series.  This  wind  tunnel,  specially  designed 
to  study  atmospheric  flow  phenomena,  incorporates  special  features 
such  as  adjustable  ceiling,  rotating  turntables,  transparent  boundary 
walls,  and  a  long  test  section  to  permit  reproduction  of  micrometeoro- 
logical  behavior  at  larger  scales.  Mean  wind  speeds  of  0.15  to  12  m/s 
can  be  obtained  in  the  EWT.  Boundary  layer  depths  1  m  thick  over  the 
downstream  6  meters  can  be  obtained  with  the  use  of  the  vortex  generators 
at  the  test  section  entrance  and  surface  roughness  on  the  floor.  The 
flexible  test  section  roof  on  the  EWT  is  adjustable  in  height  to  permit 
the  longitudinal  pressure  gradient  to  be  set  at  zero.  The  vortex 
generators  at  the  tunnel's  entrance  were  followed  by  10  m  of  smooth 
floor,  and  a  3  m  approach  ramp  to  the  1:85  scaled  topography  at  the 
China  Lake  site. 

3.2  Model 

A  1:170  scale  model  of  the  China  Lake  topography  was  constructed 
for  use  in  the  Meteorological  Wind  Tunnel  and  a  1:85  scale  model  was 
constructed  for  use  in  the  Environmental  Wind  Tunnel.  The  topographic 
relief  of  the  China  Lake  site  is  shown  in  Figure  6.  Both  these  models 
were  constructed  of  0.64  cm  thick  styrofoam  sheets.  A  cylindrical 
plenum  manufactured  with  perforated  upper  plate  was  centered  in  the 
middle  of  the  test  site  pond.  The  source  gas.  Argon,  stored  in  a  high 
pressure  cylinder  was  directed  through  a  solenoid  valve,  a  flow  meter, 
and  onto  the  circular  area  source  mounted  in  the  model  pond.  Typical 
curves  of  the  boiloff  duration  and  source  diameter  are  shown  in  Figure  7. 
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All  source  release  conditons  were  step  functions;  thus,  their  profiles 
can  be  recreated  from  the  data  in  Table  3. 

3.3  Flow  Visualization  Techniques 

Smoke  was  used  to  define  plume  behavior  over  the  China  Lake  site. 

The  smoke  was  produced  by  passing  the  simulation  gas,  Argon,  through  a 
container  of  titanium  tetrachloride  located  outside  the  wind  tunnel. 

The  plume  was  illuminated  with  arc- lamp  beams.  A  visible  record  was 
obtained  by  means  of  pictures  taken  with  a  Speed  Graphic  camera  utilizing 
Polaroid  film  for  immediate  examination.  Additional  still  pictures 
were  obtained  with  a  35  mm  camera.  The  color  motion  pictures  were  taken 
with  a  Bolex  motion  picture  camera. 

3.4  Wind  Profiles  and  Turbulence  Measurements 

Velocity  profile  measurements  and  reference  wind  speed  conditions 
were  obtained  with  a  Thermo-Systems  Inc.  (TSI)  1050  anemometer  and  a 
TSI  model  1210  hot  film  probe.  Turbulence  measurements  were  made  with 
this  system  for  the  longitudinal  velocity  component  and  with  a  TSI  split 
film  probe  connected  to  two  TSI  1CS0  anemometers  for  both  longitudinal 
and  vertical  component  measurements.  Since  the  voltage  response  of 
these  anemometers  is  non-linear  uith  respect  to  velocity,  a  multi-point 
calibration  of  system  response  v.rsus  velocity  was  utilized  for  data 

i 

reduction. 

The  velocity  standard  util:'  :ed  in  the  present  study  was  that 
depicted  in  Figure  8.  This  comisted  of  a  Matheson  model  8116-0154 
mass  flowmeter,  a  Yellowsprings  thermistor,  and  a  profile  conditioning 
section  constructed  by  the  Engii  Bering  Research  Center  shop.  The  mass 
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flowmeter  measures  mass  flow  rate  independent  of  temperature  and 
pressure,  the  thermistor  measures  the  temperature  at  the  exit  conditions, 
and  the  profile  conditioning  section  forms  a  flat  velocity  profile 
of  very  low  turbulence  at  the  position  where  the  probe  is  to  be  located. 
Incorporating  a  measurement  of  the  ambient  atmospheric  pressure  and  a 
profile  correction  factor  permits  the  calibration  of  velocity  at  the 
measurement  station  from  0. 0-2.0  m/s  + 5.0  cm/s. 

During  calibration  of  the  single  film  anemometer,  the  anemometer 
voltage  response  values  over  the  velocity  range  of  interest  were  fit  to 
an  expression  similar  to  that  of  King's  law  (Sandbom,  1972)  but  with  a 
variable  exponent  determined  by  least  squares  method.  The  accuracy  of 
this  technique  is  approximately  +2  percent  of  the  actual  longitudinal 
velocity. 

The  split  film  probe  was  mounted  on  a  rotatable  mechanism  with  a 
precision  protractor  accurate  to  + 2  minutes  and  positioned  in  front  of 
the  velocity  standard.  The  calibration  procedure  is  described  in  the 
TSI  Technical  Bulletin  20.  The  suggested  calibration  equations  were 
fit  to  the  data  by  a  least  squares  method.  The  accuracy  of  this  technique 
is  approximately  +5  percent  of  the  actual  longitudinal  or  vertical 
velocity  values. 

The  velocity  sensors  were  mounted  on  a  vertical  traverse  and  positioned 
over  the  measurement  location  on  the  model.  The  anemometer's  responses 
were  fed  to  a  Preston  analog-to-digital  converter  and  then  directly  to 
a  HP-1000  minicomputer  for  immediate  interpretation.  The  HP-1000 
computer  also  controls  probe  position.  A  flow  chart  depicting  the  control 
sequence  for  this  process  is  presented  in  Figure  9. 
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3.5  Concentration  Measurements 

The  concentrations  of  methane  produced  during  an  LNG  spill  are 
inherently  time  dependent.  It  is  necessary  to  have  a  frequency  response 
to  concentration  fluctuations  of  at  least  50  Hz  to  isolate  peaks  of 
methane  concentrations  above  5  percent  (the  lower  flammability  limit 
of  methane  in  air,  LFL) ;  hence,  an  aspirating  hot  film  probe  was 
used  for  this  study. 

3.5.1  Hot  Film  Aspirating  Probe 

The  basic  principles  governing  the  behavior  of  such  a  probe  have 
been  discussed  by  Blackshear  and  Fingerson  (1962),  Brown  and  Rebollo 
(1972),  and  Kuretsky  (1967).  A  diagram  of  the  design  of  this  probe  is 
presented  in  Figure  11.  A  vacuum  source  sufficient  to  choke  the  flow 
through  the  small  orifice  just  downwind  of  the  sensing  elements  was 
applied.  Only  one  of  the  two  films  in  this  special  probe  was  an  active 
element  for  the  measurement  of  concentration  in  the  present  study.  This 
film  was  operated  in  a  constant  temperature  mode  at  a  temperature  above 
that  of  the  ambient  air  temperature.  A  feedback  amplifier  maintained 
a  constant  overheat  resistance  through  adjustment  of  the  heating  current. 

A  change  in  output  voltage  from  this  sensor  circuit  corresponds  to  a 
change  in  heat  transfer  between  the  hot-wire  and  the  sampling  environment. 

The  heat  transfer  rate  from  a  hot  cylindrical  film  to  a  gas  flowing 
over  it  depends  primarily  upon  the  film  diameter,  the  temperature  difference 
between  the  film  and  the  gas,  the  thermal  conductivity  and  viscosity  of 
the  gas,  and  the  gas  velocity.  For  a  film  in  an  aspirated  probe  with  a 
sonic  throat,  the  gas  velocity  can  be  expressed  as  a  function  of  the  ratio 


of  the  probe  cross-sectional  area  at  the  film  position  to  the  area  at 
the  throat,  the  specific  heat  ratio,  and  the  speed  of  sound  in  the 
gas.  The  latter  two  parameters,  as  well  as  the  thermal  conductivity  and 
viscosity  of  the  gas  mentioned  earlier,  are  determined  by  the  gas 
composition  and  temperature .  Hence,  for  a  fixed  probe  geometry  and 
film  temperature,  the  heat  transfer  rate,  or  the  related  voltage  drop 
across  the  film  is  a  function  of  only  the  gas  composition  and  temperature. 
Since  all  tests  performed  in  this  study  were  in  an  isothermal  flow 
situation  the  film's  response  was  only  a  function  of  gas  composition. 

During  probe  calibration  known  compositions  of  Argon-air  mixtures 
were  passed  through  a  pre-heat  exchanger  to  condition  the  gas  to  the 
tunnel  temperature  environment.  These  known  compositions  were  produced 
from  a  bottle  of  pure  Argon  and  bottle  of  pure  air  passed  through  a 
Matheson  gas  proportioner  or  drawn  from  a  bottle  of  prepared  gas 
conposition  provided  by  Matheson  Laboratories.  Figure  10  displays 
the  measured  variation  of  the  voltage  drop  with  percentage  of  Argon 
in  an  Argon-air  mixture  fc c  three  different  values  of  the  film  temper¬ 
ature  overheat.  For  an  overheat  ratio  (temperature  of  film/ ambient 
temperature)  of  1.75  the  voltage  drop  varies  linearly  with  Argon  concen¬ 
tration  and  has  the  maximum  sensitivity.  This  particular  overheat  ratio 
was  used  during  all  wind  tunnel  measurements. 

3.5.2  Errors  in  Concentration  Measurement 

The  effective  sampling  area  of  the  probe  inlet  is  a  function  of 
the  probe's  aspiration  rate  and  the  distribution  of  approach  velocities 
of  the  gases  to  be  sampled.  A  calculation  of  the  effective  sampling 
area  during  all  tests  suggests  that  the  effective  sampling  area  was 
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2 

always  less  than  the  area  of  the  probe's  inlet,  1.88  cm  .  Thus  the 

resolution  of  the  concentration  measurements  as  applied  to  the  China 

2 

Lake  site  is  *1.6  m  . 

The  travel  time  from  the  sensor  to  the  sonic  choke  limits  the 
upper  frequency  response  of  the  probe.  At  high  frequencies  the 
correlation  between  concentration  fluctuation  and  velocity  fluctuations 
(velocity  fluctuations  are  a  result  of  the  changes  of  sonic  velocity 
with  concentration)  at  the  sensor  begin  to  decline.  The  CSU  aspirated 
probe  is  expected  to  have  a  1000  Hz  upper  frequency  response,  but,  to 
inf>rove  signal  to  noise  characteristics,  the  signal  was  filtered  at 
200  Hz.  This  is  well  above  the  frequencies  of  concentration  fluctua¬ 
tions  that  were  expected  to  occur. 

The  errors  caused  by  a  linearity  assumption  in  the  reduction  of 
concentration  data  are  approximately  the  component  value  (percent 
Argon)  +0. 75  percent.  The  errors  caused  by  calibration  change  due  to 
tenqjerature  drift  are  approximately  0.1  percent  of  the  component  value 
per  degree  centigrade.  Since  the  tunnel  temperatures  vary  at  most 
♦5°C  during  a  given  test  period  the  maximum  error  due  to  temperature 
drift  would  be  0.5  percent  of  the  component  value. 


4.0  TEST  PROGRAM  RESULTS 

The  test  program  consisted  of  four  different  test  series.  The 

first  test  series  objectives  were: 

•  To  determine  the  general  behavior  of  LNG  spills  at  the  China 
Lake  facility. 

•  To  obtain  concentrations  and  photographs  of  the  behavior  of  these 
tests  scaled  at  1:170  to  compare  with  tests  performed  at  a  scale 
of  1:85  as  a  measure  of  the  quality  of  the  simulation  criteria. 

•  To  examine  the  scaling  implications  of  adjusting  the  initial 
specific  gravity  of  the  model  to  be  different  from  that  of  the 
prototype. 

The  second  test  series  objectives  were: 

•  To  determine  the  general  behavior  of  LNG  spills  at  the  China 
Lake  facility. 

•  To  obtain  concentrations  and  photographs  of  the  behavior  of  these 
tests  scaled  at  1:85  to  compare  to  the  1:170  scales  test 

.  performed  earlier. 

The  third  test  series  objectives  were: 

•  To  obtain  concentration  data  for  the  field  tests,  LNG-18  and 
LNG-19  at  a  scale  of  1:85  in  the  presence  of  the  China  Lake 
topography  and  with  a  smooth  floor. 

•  To  obtain  spectral  measurements  of  longitudinal  and  vertical 
velocity  components  at  two  different  sites  on  the  China  Lake  model. 

The  fourth  test  series  objectives  were: 

•  To  simulate  the  field  tests,  LNG-18,  19,  20,  and  21  and  obtain 
concentration  measurement  of  these  tests  over  the  1:85  scale  model 


of  China  Lake. 
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A  summary  of  all  tests  simulated  in  the  laboratory  is  presented  to 
Table  1.  A  detailed  summary  of  the  prototype  and  model  test  conditions 
for  LNG-18,  19,  20,  and  21  is  presented  in  Tables  2  and  3  respectively. 

All  dimensions  reported  in  the  following  sections  have  been  converted 
to  equivalent  full-scale  values  appropriate  to  the  China  Lake  site. 

The  origin  is  referenced  as  the  LNG  spill  point.  The  positive  x-axis 

is  in  the  direction  of  the  prevailing  wind  for  all  coordinate  systems  unless 

otherwise  mentioned. 

4.1  Characteristics  of  the  Approach  Velocities 

Measurements  of  the  approach  flow  characteristics  were  obtained  for 
the  modeled  flow  over  the  China  Lake  scale  topography.  As  discussed 
in  Section  2.1.1  these  characteristic  length  and  velocity  scales  should 
be  comparable  to  those  expected  to  occur  over  the  China  Lake  site. 

Counihan  (1975)  has  summarized  the  values  of  aerodynamic  roughness,  zq, 
longitudinal  velocity  integral  length  scale,  Ax»  and  the  power  law 
index,  1/n,  that  may  be  expected  to  occur  in  the  atmosphere.  Table  6 
compares  values  of  these  quantities  as  cited  by  Counihan  and  values  scaled 
up  from  the  model  tests.  Figures  13  and  14  show  the  profiles  of  mean 
velocity  and  local  turbulent  intensity  respectively. 

4.2  First  Test  Series  Results 

A  1:170  scale  model  of  the  China  Lake  site  was  placed  in  the 
Meteorological  Kind  Tunnel  facility  (MWT) .  Five  different  tests  were 
performed  on  the  model.  Each  was  performed  with  a  spill  volume  of  5.95  cubic 
meters  of  LNG  and  a  boiloff  rate  of  33.6  kg/s.  Two  different  wind 
directions  and  three  different  wind  speeds  were  considered.  These  conditions 


are  summarized  in  Table  3.  Black-and-white  Polaroid  photographs  were 
taken  of  the  different  spill  conditions  so  that  plume  configurations 
could  be  compared  to  tests  performed  over  the  1:85  scale  model.  Time 
dependent  concentration  measurements  were  obtained  at  three  axial  distances 
downwind  for  each  test.  Table  5-1  summarizes  the  times  of  arrival,  peak 
concentrations,  time  for  passage  of  the  plume  and  value  of  peak  concentra¬ 
tion  for  the  various  measurements  sites. 

Visualization  of  the  lowest  wind  speed  simulated,  2.2  m/s, 
indicated  an  unexpectedly  large  vertical  dispersion  rate.  The  concen¬ 
tration  measurements  confirmed  this  observation.  It  was  concluded 
that  this  abnormal  vertical  variance  was  associated  with  instabilities 
in  the  MWT  while  operating  below  a  minimum  stable  operational  wind 
speed. 

Tests  were  also  performed  to  test  the  advisability  of  relaxing 
precise  equality  of  the  density,  mass,  and  momentum  ratios  discussed 
in  the  simulation  chapter  (see  Section  2.2  and  2.3).  Argon  gas 
(specific  gravity  of  1.38)  ,  carbon  dioxide  gas  (specific  gravity  of 
1.52),  and  compressed  air  (specific  gravity  of  1.00)  were  used  as 
source  gases  released  from  an  area  source  on  a  smooth  floor.  Wind 
speeds  were  low  enough  so  that  pronounced  gravity  spreading  was 
exhibited.  From  these  tests  in  which  both  visual  plume  pattern  and 
concentration  measurements  were  obtained  it  was  concluded  that  very 
similar  behavior  was  exhibited  by  both  simulation  schemes.  A  +d5% 
variation  in  source  Froude  number  was  not  detectable  in  the  resultant 


4.3  Second  Test  Series  Results 

A  1:85  scale  model  of  the  China  Lake  site  was  placed  in  the  Environ¬ 
mental  Wind  Tunnel  facility  (EWT) .  Four  different  tests  were  performed 
on  the  model.  The  transient  boiloff  experiments  were  adjusted  to  simu¬ 
late  a  spill  volume  of  5.95  cubic  meters  of  LNG.  Steady  state  boiloff 
experiments  were  run  at  a  continuous  boiloff  rate  of  33.6  kg/s.  The 
time  duration  of  the  boiloff  was  for  practical  purposes  infinite.  Three 
different  wind  speeds  were  tested  approaching  from  the  southwest. 

These  conditions  are  summarized  in  Table  3.  Black-and-white  Polaroid 
photographs  and  color  slides  were  obtained  for  the  steady  boiloff 
releases  and  a  16  mm  color  movie  was  prepared  for  all  tests.  Table  5-2 
and  5-3  summarize  the  time  of  arrival,  peak  concentrations,  time  for 
passage  of  the  plume  and  the  value  of  peak  concentrations. 

Measurements  from  the  first  test  series  in  the  MWT  on  a  1:170  scale 
model  and  measurements  from  this  test  series  on  a  1:85  scale  model  may 
be  compared.  Inspecting  Table  4,  it  can  be  concluded  that  concentrations 
measured  at  different  model  scales  are  generally  of  the  same  order  except 
for  the  cases  in  which  the  MWT  tunnel  was  run  at  its  lowes"  speed.  It 
was  mentioned  in  Section  4.2  that  this  speed  was  felt  to  be  below  the  opera¬ 
tional  limit  of  the  MWT.  Since  the  number  of  measurement  locations  and  test 
replications  were  very  limited  in  the  MWT  no  conclusions  can  be  made  con¬ 
cerning  the  quantitative  influence  of  model  scale  on  the  simulation. 

4.4  Third  Test  Series  Results 

The  China  Lake  boiloff  rate,  boiloff  duration,  and  wind  speed  of  the 
LNG-18  and  LNG-19  tests  were  simulated  in  the  Environmental  Wind  Tunnel 
with  a  smooth  floor.  The  same  tests  were  then  repeated,  but  this  time 


the  topography  of  the  China  Lake  site  was  included!  Thorough  concen¬ 
tration  measurements  downwind  of  these  four  tests  were  obtained.  A 
summary  of  the  test  conditions  for  these  four  tests  is  presented  in 
Table  3.  Approximate  hand  drawn  ground  level  peak  concentration  contours 
are  shown  in  Figure.'  15-1  through  15-4  for  all  four  te.  ts.  The  maximum 
normalized  concentration  'oefficient  versus  downwind  distances  for  the 
two  tests  without  the  topo/raphy  is  presented  in  Figure  16-1.  A  summary 
of  the  times  of  arrival,  pe  ik  concer  ration  and  passage  of  the  plume 
and  the  maximum  peak  concentration  observed  is  presented  in  Table  5-4. 

It  is  see.t  from  a  comparison  uf  the  similar  tests  one  performed  with 
topography  and  one  performed  on  a  smooth  floor,  that  t  le  dispersion  at 
the  China  Lake  site  is  greater  than  what  would  be  expe *ienced  if  the 
spill  occurred  in  a  very  smooth  and  flat  area. 

Spectral  measurements  of  the  longitudinal  and  vertical  velocity 
fluctuations  were  obtained  at  two  different  locations  at  a  height  of 
2  m  over  the  China  Lake  scaled  topography.  A  discussion  of  these  measure¬ 
ments  appears  in  Section  4.1  on  the  simulated  characteristics  of  wind 
motion  over  the  China  Lake  site. 

4.5  Fourth  Test  Series  Results 

A  i :85  scale  model  of  the  China  Lake  site  was  placed  in  the 
Environmental  Wind  Tunnel.  Concentration  measurements  were  obtained 
downwind  of  the  simulated  field  tests  LNG-18,  LNG-19,  LNG-20,  and  LNG-21. 

A  summary  of  the  field  conditions  that  were  simulated  is  presented  in 
Table  1.  A  summary  of  the  model  conditions  for  these  tests  is  presented 
T 


Unfortunately  the  wind  directions  provided  by  the  field  investigators 
were  in  error.  These  two  tests  were  rerun  in  the  fourth  test  series. 
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in  Table  2.  Ground  level  peak  concentration  contours  are  shown  in 
Figures  15-5  through  15-8  for  each  test.  The  maximum  normalized 
concentration  versus  downwind  distance  for  all  four  tests  are  shown 
in  Figure  16-2.  A  summary  of  the  time  of  arrival,  peak  concentration 
and  passage  of  the  plume  and  the  maximum  peak  concentration  observed 
is  presented  in  Table  5-5. 

4.6  Comparison  with  Field  Data 

As  part  of  the  current  China  Lake  test  series,  field  concentration 
measurements  were  obtained  over  two  independent  measurement  grids.  The 
Naval  Weapons  Test  Center  established  a  grid  of  ten  different  concen¬ 
tration  measurement  stations  and  the  Lawrence  Livermore  Laboratory  (LLL) 
provided  eight  towers  with  a  variety  of  concentration  sampling  equipment. 
The  primary  purpose  of  the  LLL  grid  was  sensor  evaluation.  Both  these 
grids  are  described  on  Figure  7. 

The  degree  to  which  physically  modeled  data  correlates  to  values 
which  are  obtained  in  the  field  is  dependent  upon  the  approximations 
which  were  assumed  in  the  formulation  of  the  model  and  the  inherent 
randomness  of  the  atmospheric  diffusion  processes.  The  assumptions 
employed  in  the  construction  of  a  physical  model  of  LNG  vapor  dispersion 
at  the  Naval  Weapons  Test  Center  are  discussed  in  Section  2.3.  The 
randomness  of  wind  directions  and  velocities  in  the  atmosphere  are  such 
that  a  single  time  realization  of  a  fixed  point  in  space  is  insuffi¬ 
cient  to  describe  the  coi.'plete  probability  distribution  of  peak  concen¬ 
trations  that  may  be  observed  at  that  point.  Without  ensemble  averaging 
of  similar  tests  in  the  field  the  values  found  during  a  single  realization 
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may  range  over  a  better  portion  of  an  unknown  probability  distribution. 
Pasquill  (1974)  notes  that  in  many  circumstances  of  practical  interest 
the  uncertainties  found  between  continuous  releases  of  gaseous  plumes 
may  at  best  be  ten  to  fifty  percent  in  the  average  and  factors  of  two 
or  more  individually.  In  addition  to  the  small  scale  effects  of  local 
randomness,  the  atmosphere  has  large  scale  effects  which  lead  to 
meandering  of  plume  mean  motion.  These  large  scale  meanderings  are 
not  modeled  in  wind  tunnels. 

The  Naval  Weapons  Test  Center  grid  consisted  of  ten  different 
concentration  sensors.  These  instruments  were  all  of  the  catalytic 
combustion  type.  The  principle  of  operation  of  these  instruments  is  that 
a  hot  catalytic  filament  causes  the  methane  to  oxidize,  and  the  rise  in 
temperature  due  to  the  reaction  changes  the  electrical  resistance  of 
the  filament.  These  detectors  are  accurate  for  only  low  (below  7%), 
slowly  varying  methane  concentration. 

Table  7  compares  peak  concentrations  observed  in  the  field  at  the 
Naval  Weapons  test  grid  points  with  those  obtained  over  a  physical 
model.  This  comparison  is  in  general  quite  poor.  There  are  several 
factors  which  may  account  for  this  scatter  in  comparable  data  over 
several  orders  of  magnitude.  That  are:  (1)  the  mean  wind  direction 
specified  for  each  wind-tunnel  test  may  have  been  in  error:  (2)  the 
fluctuations  present  during  the  field  tests  were  as  large  as  +50° 
(physical  modeling  of  large  wind  direction  fluctuations  is  not 
possible  in  a  wind  tunnel);  (3)  the  wind  speed  observed  in  the  field 
changed  by  as  much  as  +1.8  m/s  during  the  tests  (this  amount  of  fluc¬ 
tuation  can  account  for  approximately  +50%  variation  in  concentration 


values);  (4)  the  peak  concentration  fluctuations  in  the  field 
tests  were  too  rapid  for  the  catalytic  sensor  to  respond; 

(5)  the  concentrations  were  too  large  for  the  catalytic  sensors 
to  respond;  (6)  the  approximations  used  in  simulating  the  LNG 
field  test  series  were  too  weak  to  achieve  proper  simulation. 

The  Lawrence  Livermore  Laboratory  obtained  concentration  tj.me 
histories  at  a  variety  of  different  heights  on  their  eight  towers  equipped 
with  concentration  sensors.  Several  different  types  of  sensors  were 
employed.  Each  of  these  detector  responses  was  verified  by  simultaneous 
grab  bag  sampling  of  the  gases  flowing  over  the  sensor.  This  technique 
provides  an  accurate  method  of  verifying  the  different  sensors' 
response  were  correct.  The  peak  field  concentration  obtained  from  the 
lowest  sensor  elevation  at  each  of  the  eight  towers  is  summarized 
along  with  approximate  model  values  in  Table  7.  Since  concentrations 
over  the  model  were  not  obtained  at  the  Lawrence  Livermore  grid  sites, 
the  values  noted  for  model  equivalents  are  only  approximate.  These 
values  were  obtained  by  interpolation  of  the  hand  drawn  ground  level 
peak  concentration  contours  in  Figures  15-15  through  15-8.  On  these 
figures  the  circled  numbers  are  the  peak  concentrations  observed 
in  the  field  on  the  Lawrence  Livermore  grid.  Figure  17  displays  the 
time  variation  of  LNG  vapor  concentration  for  both  modeled  and  field 
data  at  LLL  grid  point  5  for  the  test  spill  LNG-21. 


The  correlation  between  Lawrence  Livermore  data  and  that  of  the 
model  data  is  generally  superior  to  that  found  between  the  Naval  Weapons 
Test  Center  data  and  that  of  the  model.  There  remain,  however,  a  number 
of  sampling  points  where  poor  agreement  exists.  Considering  each 
test  point  at  a  time  for  LNG-18  shows  reasonable  comparison,  within 
a  50%  of  the  field  value,  for  the  near  field  grid  points  1,  2,  and  3, 
and  poor  comparison  on  grid  point  four.  The  reason  for  these  discon¬ 
tinuities  in  field-model  comparisons  may  be  any  combination  of  the 
factors  mentioned  previously.  In  this  case  the  differences  appear 
to  be  caused  by  the  small  number  of  measurement  locations  in  both 
field  and  model  tests  and  the  variability  of  wind  direction  in  the 
field.  For  LNG-19  the  quality  of  comparison  between  model  and  field 
data  is  somewhat  poor.  The  decay  of  concentration  with  distance  from 
the  source  appears  to  agree,  but  the  direction  of  the  plume  appears 
to  be  different.  This  result  suggests  a  change  in  wind  direction 
between  field  data  and  what  was  modeled.  Here  again  as  in  LNG-18  an 
insufficient  number  of  model  or  field  measurement  locations  were  taken 
to  define  the  concentration  field  properly.  In  LNG-20  the  comparison 
between  field  and  model  again  appears  to  be  poor.  The  laboratory 
model  predicts  that  at  this  higher  wind  speed  (12.4  m/s  at  2  m)  the 
LNG  plume  has  very  little  lateral  spread,  but  the  field  measurement 
show  concentrations  at  large  distances  from  the  plume's  mean  axis. 

This  suggests  a  large  variation  in  wind  direction  or  an  error  in 
the  mean  wind  direction.  For  this  test  and  the  LNG-21  a 
sufficient  number  of  measurement  locations  were  used  to  define 
the  model  ground  level  contours  properly,  47  and  91  points 
respectively.  LNG-21  shows  the  best  comparison  between  model  and 
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field  results  for  the  four  tests  modeled.  Each  measurement  location 
has  an  acceptable  comparison  considering  the  field  variation  of  wind 
direction  and  velocity  and  that  there  are  insufficient  field  data 
points  to  define  ground  level  peak  concentration  patterns. 
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TABLE  1.  SUMMARY  OF  CONCENTRATION  TESTS 

(All  results  at  equivalent  field 
condition  levels) 
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TABLE  2 .  PROTO! 


CHARACTERISTICS 
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Specific  Gravity  •  boiloff  temp. 

Rind  Speed  .  U  #  2  *  (!) 

Wind  Direction 

Stability  (Pasqui 11 -Gifford  Category) 

Humidity  (%) 

*  • 

Density  Retio,  ^  ^ 

*ir  •  boiloff  temp. 

Reynolds  No.. ,  2, 

tl2 

Froude  No  J — - - \ 

\l  7  °/  ,  2  a 


Richardson  No., 


(T|0.  -  T2.)»W 

TS.(C10,  ‘  02.)2 


ror  source  release  conditi 


At  the  10  meter  tower 


uns  jer 


■  •  6“  •  c 


T.  •  lll.6J°K,  0 
b .  o . 


LNG 
-5  _2 


•  422.63  -J  .  oNC  R  boilof 


•  i.526  *  10  ■  Vs, 


1.197  it 


0m  • 


CONDITIONS 


l.N'G  Test  \o 


18 

19 

20 

21 

20 

:o 

20 

20 

4.39 

5.2 

4.S 

4  2 

1000 

1184,5 

1025 

9S6.T 

67 

S9 

-7 

S3 

27V 

37.4 

24.8 

33.6 

14.9 

20. 1 

13.3 

18.0 

1.5S 

1.55 

l .  SS 

I  .SS 

6.7 

5.1 

12.4 

4.9 

214° 

260° 

256° 

224° 

C 

C-D 

D 

c 

16 

29 

15 

21 

0.55 

0.55 

0.5S 

0.55 

I.SxlO6 

6.7»I06 

1.6xl07 

6.4X106 

0.42 

0.24 

1.42 

0.22 

-0.  IS 


86 


3 

■ 


-0.06 


0.11 


TABLE  3.  MODEL  CONDITIONS 
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TABLE  4.  COMPARISON  OF  SIMILAR  LNG  TESTS  AT  DIFFERENT  MODEL  SCALES 


X** 

Cm) 

WIND 

DIR. 

WIND 

SPEED 

(m/s) 

* 

X 

1:170 

WIND 

DIR. 

WIND 

SPEED 

(m/s) 

X1 : 85 

23.0 

SW 

2.2 

SW 

2.3 

28.0 

If 

ft 

14.7 

ft 

If 

32.5 

45.0 

11 

II 

12.5 

ft 

ft 

21.5 

63.0 

If 

tl 

5.2 

tt 

ft 

68.5 

tl 

It 

ft 

It 

11.3 

91.5 

II 

II 

ft 

ft 

10.4 

23.0 

SW 

8.9 

SW 
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28.0 

If 

tl 
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It 

tl 
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It 

tl 

11.3 

II 

It 

20.1 

63.0 

ft 

If 

8.0 

tt 

tt 

68.5 

If 

II 

It 

ft 

8.4 

91.5 

If 

II 

It 

tt 

7.0 

23.0 

w 

2.2 

SW 

2.3 

28.0 

tl 

tt 

23.1 

ft 

tt 

32.5 

45.0 

ft 

ft 

14.7 

ft 

tt 

21.5 

63.0 

ft 

It 

12.5 

tt 

ft 

63.0 

It 

II 

tl 

It 

11.3 

91.5 

tl 

II 

tl 

If 

10.4 

‘Data  taken  with  hot  film  aspirating  probe 

** (y  -  o,  z  =  o) 
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TABLE  5-1.  SUMMARY  OF  CONCENTRATION  DATA 


Position 

x  y  z 

(m)  (m)  (m) 

Maximum 

Peak 

Concentration 

(%) 

Time  of 
Arrival 
(sec) 

Time  of 
Peak 
(sec) 

Approximate 
Time  of 
Passage 
(sec) 

Test -Run 

No.  1-1 

0 

14.7 

36.9 

147.6 

393.6 

12.5 

49.2 

159.9 

528.9 

0 

5.2 

61.5 

203.0 

639.6 

Test-Run 

No.  1-2 

27.0 

26.0 

24.6 

116.9 

184.5 

0 

14.7 

18.5 

86.1 

258.3 

61.0 

0 

13.6 

24.6 

129.2 

393.6 

Test -Run 

No.  1-3 

0 

23.1 

18.5 

98.4 

147.6 

44.0 

0 

11.3 

18.5 

86.1 

147.6 

61.0 

0 

8.0 

12.3 

73.8 

110.7 

Test-Run 

No.  1-4 

27.0 

0 

0 

23.1 

12.3 

49.2 

221.4 

0 

14.7 

12.3 

49.2 

215.3 

61.0 

0 

0 

12.5 

43.1 

159.9 

282.9 

Test -Run 

No.  1-5 

0 

26.0 

24.6 

98.4 

196.8 

0 

0 

19.0 

30.8 

110.7 

209.1 

13.6 

36.9 

110.7 

178.4 
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TABLE  5-2.  SUMMARY  OF  CONCENTRATION  DATA 


* 

X 

(m) 

Position 

y 

(nO 

(m) 

Maximum  Approximate 

Peak  Time  of  Time  of  Time  of 

Concentration  Arrival  Peak  Passage 

(%)  (sec)  (sec)  (sec) 

Test-Run 

No.  2-1  (Wind  Direction 

o  * 

=  225°) 

29.9 

0 

0 

32.5 

26.7 

91.5 

213.5 

54.9 

-17.1 

0 

11.8 

61.0 

160.2 

266.9 

45.8 

0 

0 

21.5 

30.5 

110.6 

278.4 

68.9 

0 

0 

11.3 

22.9 

129.6 

373.7 

91.5 

0 

0 

10.4 

95.1 

152.5 

274.5 

42.1 

-17.7 

0 

19.4 

30.5 

99.1 

266.9 

63.4 

-26.2 

0 

12.0 

76.3 

125.9 

240.3 

42.1 

17.7 

0 

26.0 

45.8 

133.5 

213.5 

j  Test-Run 

No.  2-2  fWind 

Direction 

»  225°) 

29.9 

0 

0 

43.5 

54.9 

-17.1 

0 

22.3 

45.8 

0 

0 

28.2 

68.9 

0 

0 

19.2 

91.5 

0 

0 

14.7 

42.1 

-17.7 

0 

28.4 

63.4 

-26.2 

0 

16.7 

Test-Run 

No.  2-5  (Wind 

Direction 

o  * 

-  225  ) 

29.9 

0 

0 

24.7 

7.6 

80.1 

137.3 

54.9 

-17.1 

0 

0.5 

72.5 

80.1 

106.8 

45.8 

0 

0 

20.1 

11.5 

61.0 

137.3 

68.9 

0 

0 

8.4 

26.7 

64.9 

114.4 

91.5 

0 

0 

7.0 

45.8 

68.6 

129.6 

42.1 

-17.7 

0 

0.3 

72.5 

76.3 

80.1 

42.1 

17.7 

0 

21.3 

22.9 

53.4 

129.6 

63.4 

26.2 

0 

14.3 

0 

0 

0 

Test-Run 

No.  2-6  (Wind 

Direction 

o  * 

-  225  ) 

29.9 

0 

27.1 

54.9 

-17.1 

0 

1.1 

45.8 

0 

21.3 

68.9 

0 

12.0 

91.5 

0 

0 

6.7 

42.1 

-17.7 

0 

1.6 

-26.2 

0 

0.5 

42.1 

17.7 

0 

24.1 

♦Coordinate  system  is  right  handed  and  the  positive  x-axis  is 
orientated  in  the  direction  of  the  wind.  Wind  direction  is  reference 
to  direction  wind  comes  from. 


Position 


x 


(■) 


3 


3 


y 

(■) 


(■) 


Maximum 

Peak 

Concentrat ion 

(*) 


Time  of  Time  of 
Arrival  Peak 

(sec)  (sec) 


0  0 

8.S  0 

8.5  0 

0  0 

5.5  0 

0  0 

8.5  0 

5.5  0 

0  0 

0  0 

8.5  0 

8.5  0 

0  0 

17.0  0 

17.0  0 

0  0 

5.5  0 

8.5  0 

0  0 

0  0 

1.5  0 

3.5  0 

0  0 

0  0 

1.5  0 

8.5  0 

0  0 

r.o  0 

1.5  0 

r.O  0 

0  0 

1.5  0 

1.5  0 

0  0 

0  0 

1.5  0 

1.5  0 

0  0 


53.9 
49.0 

33.9 
29.0 
31.2 

27.1 

27.2 

32.5 
25.  S 

22.6 

17.0 

17.1 

20.5 

12.0 

0.46 

15.3 

15.1 

17.8 

15.1 

12.4 

11.7 

15. 7 
11.0 
12.  S 

9.8 
8.6 
8.2 
8.0 
1.6 
2.1 

10.9 
9.0 

9.5 

8.9 
9.3 

6.9 

7.7 

5.9 


9.6 

65.0 

9.6 

68.9 

13.4 

68.9 

9.6 

51.6 

11.5 

44.0 

13.4 

51.6 

13.4 

76.5 

15.3 

74.6 

17.2 

66.9 

17.2 

65.0 

13.4 

68.9 

17.2 

66.9 

15.3 

61.2 

28.1 

49.1 

47.8 

47.8 

19.1 

63.1 

17.2 

65.0 

19.1 

76. S 

19.1 

72.7 

19.1 

76.5 

19.1 

84.2 

21.0 

72.7 

msem 

24.9 

26.8 

57.4 

28.7 

93.7 

30.6 

88.0 

49.7 

53.6 

44.0 

23.0 

88.0 

28.7 

86.1 

42.1 

80.3 

30.6 

91.8 

36.3 

51.6 

78.4 

30.6 

72.7 

34.4 

59.3 

Approximate 
Time  of 
Passage 

(sec) 


162.6 

99. 5 

95.6 
65. 0 

114.8 

166.4 

114.8 

139.6 
153.0 
172.1 

118.6 

122.4 

147.3 
137.7 

47.8 

172.1 

137.7 
1 

1  /  ** .  i 

153.0 

120. 5 

153.0 

133.9 
172.1 

124.3 

114.8 

122.4 

133.9 
65.0 
57.4 

154.9 

158.8 

139.6 

166.4 

135.8 

114.8 

133.9 
133.9 
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TABLE  S-3 .  (continued) 
Test -Run  No.  3-2 


Position 

Maximum 

Approximate 

Peak 

Tine  of 

Time  of 

Time  of 

X 

Y 

z 

Concentration 

Arrival 

Peak 

Passage 

(■) 

(■) 

(■) 

(%) 

(sec) 

(sec) 

(sec) 

18.  S 

0 

0 

48.7 

9.6 

67.0 

136.0 

27.0 

0 

0 

42.8 

15.3 

70.9 

105.3 

27.0 

12.8 

0 

38.  S 

11.5 

51.7 

111.1 

27.0 

25.5 

0 

0.0 

0 

0 

0 

27.0 

-12.8 

0 

9.9 

53.6 

57.5 

61.3 

27.0 

-25. S 

0 

0.0 

0 

0 

0 

44.0 

0 

0 

34.5 

21.1 

53.6 

149.4 

61.0 

0 

0 

24.4 

19.2 

72.8 

134.1 

61.0 

12.8 

0 

25.2 

19.2 

74.7 

141.7 

61.0 

25.5 

0 

l.SS 

26.8 

28.7 

30.6 

61.0 

-12.8 

0 

18.6 

63.2 

76.6 

86.2 

61.0 

-25.5 

0 

0 

0 

0 

0 

78.0 

0 

0 

22.1 

23.0 

68.9 

141.7 

9S.0 

0 

0 

18.3 

38.3 

76.6 

137.9 

112.0 

0 

0 

14.5 

32.6 

80.4 

143.6 

112.0 

12.8 

0 

13.7 

28.7 

95.8 

162.8 

112.0 

25.  S 

0 

6.8 

44.0 

80.4 

107.2 

112.0 

-12.8 

0 

6.1 

SI. 7 

84.3 

120.6 

112.0 

-25.5 

0 

0.0 

0 

0 

0 

129.0 

0 

0 

12.4 

38.3 

86.2 

153.2 
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TABLE  5-3.  (continued) 


Test-Run  No.  3-3  (Wind  Direction  =  203°)* 


Position 


Maximum 

Peak 


Approximate 

Time  of  Time  of  Time  of 


X 

y 

z 

Concentration 

(■) 

(■) 

(■) 

(%) 

22. S 

0 

0 

30.9 

22.  S 

8.5 

0 

22.5 

22.5 

17.0 

0 

5.8 

22.5 

-8.5 

0 

19.6 

22.5 

-17.0 

0 

4.9 

22.5 

-2S.S 

0 

0.2 

35 . 3 

0 

0 

13.2 

3S.3 

12.8 

0 

14.8 

43.8 

0 

0 

13.8 

43.8 

-17.0 

0 

6.7 

48.0 

0 

0 

12.9 

48.0 

8.5 

0 

15.8 

48.0 

17.0 

0 

4.6 

48.0 

25.5 

0 

1.6 

48.0 

-8.5 

0 

9.0 

48.0 

-17.0 

0 

4.6 

48.0 

-2S.5 

0 

1.0 

56.5 

4.3 

0 

8.18 

60.8 

0 

0 

6.7 

73.5 

0 

0 

7.3 

73.5 

8.5 

0 

7.3 

73. S 

17.0 

0 

4.0 

73.5 

25.5 

0 

1.6 

73.5 

34.0 

0 

0.1 

73.5 

-8.5 

0 

S.8 

73.5 

-17.0 

0 

3.4 

73.5 

-25.5 

0 

3.7 

73.5 

-34.0 

0 

0.1 

86.3 

0 

0 

6.1 

86.3 

8.5 

0 

S.8 

86.3 

17.0 

0 

4.3 

86.3 

25.5 

0 

2.7 

86.3 

34.0 

0 

1.1 

86.3 

-8.5 

0 

6.  1 

86.3 

-17.0 

0 

1.9 

107.  S 

0 

0 

7.0 

107.5 

8.S 

0 

5.1 

107.5 

17.0 

0 

6.1 

107.5 

25.5 

0 

1.9 

107.5 

34.0 

0 

1.1 

107.5 

-8.S 

0 

4.2 

107.5 

-17.0 

0 

2.2 

Arrival  Peak  Passage 

(sec)  (sec)  (sec) 


11.  S 

38.3 

141.5 

15.3 

42.1 

97.6 

45.9 

89.9 

93.7 

15.3 

44.0 

124.3 

32.5 

36.3 

42.1 

0 

0 

0 

23.0 

65.0 

118.6 

34.4 

4S.9 

120.5 

17.2 

68.9 

137.7 

30.6 

80.3 

114.8 

19.1 

59.3 

130.1 

23.0 

86.1 

133.9 

19.1 

59.3 

89.9 

66.9 

68.9 

78.4 

23.0 

68.9 

120. S 

23.0 

61.2 

93.7 

42.1 

44.0 

49.7 

19.1 

47.8 

97.6 

26.8 

57.4 

110.9 

24.9 

49.7 

114.8 

23.0 

57.4 

137.7 

24.9 

36.3 

78.4 

45.9 

66.9 

89.9 

0 

0 

0 

24.9 

80.3 

114.8 

24.9 

51.6 

86.1 

61.2 

72.7 

118.6 

0 

0 

0 

24.9 

93.7 

137.7 

23.0 

49.7 

133.9 

23.0 

99.5 

122.4 

51.6 

59.3 

105.2 

66.9 

80.3 

88.0 

23.0 

61.2 

126.2 

45.9 

80.3 

114.8 

23.0 

70.8 

154.9 

23.0 

74.6 

137.7 

44.0 

72.7 

141.5 

45.9 

95.6 

105.2 

36.3 

51.6 

57.4 

49.7 

84.2 

149.2 

51.6 

107.1 

137.7 

‘Coordinate  system  is  right  handed  and  the  positive  x-axis  is 
orientated  in  the  direction  of  the  wind.  Wind  direction  is  reference 
to  direction  wind  comes  from. 
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TABLE  5-3.  (continued) 

Test-Run  No.  3-3  (continued)  Wind  Direction  *  203°)* 


X 

(■) 

Position 

y 

(■) 

z 

(■) 

Maximus 

Peak 

Concentration 

(») 

Time  of 
Arrival 
(sec) 

Time  of 
Peak 

(sec) 

Approximate 
Time  of 
Passage 

(sec) 

107.  S 

-2S.S 

0 

1.1 

49.7 

72.7 

80.3 

120.3 

0 

0 

7.4 

30.6 

74.6 

141.  S 

120.3 

8.S 

0 

6.4 

36.3 

68.9 

145.4 

120.3 

17.0 

0 

4 .  S 

30.6 

82.2 

130.1 

120.3 

2S.S 

0 

3.0 

28.7 

101. 4 

120.5 

120.3 

34.0 

0 

1.7 

68.9 

86.1 

103.3 

120.3 

42. S 

0 

O.S 

76. S 

76.  S 

76.5 

120.3 

-8.5 

0 

S.4 

28.7 

74.6 

147.3 

120.3 

-17.0 

0 

1.7 

34.4 

76. S 

126.2 

120.3 

-25. 5 

0 

0.3 

0 

0 

0 

14S.S 

0 

0 

3.2 

28.7 

66.9 

137.7 

14S.  8 

17.0 

0 

2.S 

34.4 

86.1 

133.9 

14S.  8 

34.0 

0 

0.0 

0 

0 

0 

♦Coordinate  system  is  right  handed  and  the  positive  x-axis  is 
orientated  in  the  direction  of  the  wind.  Wind  direction  is  reference 
to  direction  wind  comes  from. 
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TABLE  5-3.  (continued) 

Test-Run  No.  3-4  (Wind  Direction  =  233°)* 


Position 


Maximus 

Peak 

Concentration 


Tine  of 
Arrival 


Tine  of 

Peak 


Approximate 
Time  of 
Passage 


(S) 

(■) 

(■) 

(%) 

(sec) 

(sec) 

(sec) 

22.5 

0 

0 

28.2 

17.2 

53.6 

113.0 

22.5 

8.5 

0 

18. 5 

13.8 

76.6 

113.0 

22.5 

17.0 

0 

2.0 

23.0 

24.9 

28.7 

22.5 

-8.S 

0 

37.6 

11.5 

57.  S 

111.1 

22.5 

-17.0 

0 

7.4 

24.9 

67.0 

74.7 

22.5 

-25.5 

0 

0.2 

0 

0 

0 

3S.3 

0 

0 

27.2 

23.0 

76.6 

115.0 

3S.3 

0 

0 

28.4 

15.3 

65.1 

145.5 

48. 0 

0 

0 

17.5 

15.3 

452.0 

153.2 

60.8 

0 

0 

17.0 

19.6 

78.5 

143.6 

60.8 

8.5 

0 

12.6 

38.3 

67.0 

113.0 

60.8 

17.0 

0 

8.2 

46.0 

88.1 

124.5 

60.8 

25. S 

0 

0.0 

0 

0 

0 

60.8 

-8.5 

0 

15.3 

21.1 

70.9 

137.9 

60.8 

-17.0 

0 

12.9 

19.6 

55.5 

105.3 

60.8 

-25.5 

0 

6.8 

26.8 

30.6 

63.2 

60.8 

-34.0 

0 

0.5 

46.0 

47.9 

47.9 

73.  S 

0 

0 

14.6 

21.1 

76.6 

118.7 

86.3 

0 

0 

13.0 

24.9 

72.8 

134.1 

86.3 

a.s 

0 

12.5 

28.7 

76.6 

172.4 

86.3 

17.0 

0 

8.7 

38.3 

84.3 

114.9 

86.3 

25.5 

0 

5.4 

42.1 

88.1 

143.6 

86.3 

34.0 

0 

0.3 

0 

0 

0 

86.3 

-8.S 

0 

10.3 

30.6 

57.5 

134.1 

86.3 

-17.0 

0 

5.2 

42.1 

63.2 

120.6 

86.3 

-2S.S 

0 

1.4 

80.4 

84.3 

103.4 

86.3 

-34.0 

0 

0.2 

0 

0 

0 

99.0 

0 

0 

9.0 

24.9 

82.3 

149.4 

99.0 

8.5 

0 

5.1 

34.5 

91.9 

157.0 

99.0 

17.0 

0 

7.8 

46.0 

61.3 

137.9 

99.0 

2S.S 

0 

4.6 

49.8 

65.1 

141.7 

99.0 

34.0 

0 

1.1 

61.3 

74.7 

80.4 

99.0 

-8.5 

0 

9.4 

26.8 

72.8 

136.0 

99.0 

-17.0 

0 

8.1 

47.9 

61.3 

134.1 

99.0 

-25.5 

0 

3.9 

80.4 

88.1 

120.6 

99.0 

-34.0 

0 

0.3 

0 

0 

0 

111.8 

0 

0 

7.6 

32.6 

72.8 

151.3 

124.  S 

0 

0 

7.9 

26.8 

76.6 

153.2 

124.  S 

8.5 

0 

7.4 

38.3 

80.4 

172.4 

124.5 

17.0 

0 

1.1 

91.9 

99.6 

126.4 

124. S 

25.5 

0 

0.8 

49.8 

99.6 

126.4 

124.  S 

-8.S 

0 

8.7 

46.0 

84.3 

153.2 

124. S 

-17.0 

0 

5.8 

34.5 

76.6 

134.1 

124.  S 

-25. S 

0 

0.3 

0 

0 

0 

♦Coordinate  system  is  right  handed  and  the  positive  x-axis  is 
orientated  in  the  direction  of  the  wind.  Wind  direction  is 
reference  to  direction  wind  comes  from. 


TABLE  5-3. (continued) 

Test -Run  No.  3-4  (continued) (Wind  Direction  *  233°) 


Position 


Maximum 

Peak 

Concentration 

(%) 


Time  of 
Arrival 
(sec) 


Time  of 
Peak 
(sec) 


Approximate 
Time  of 
Passage 
(sec) 


137.3 

137.3 

137.3 

137.3 

137.3 

137.3 

137.3 

137.3 


♦Coordinate  system  is  right  handed  and  the  positive  x-axis  is  orientated 
in  the  direction  of  the  wind.  Wind  direction  is  reference  to  direction 
wind  comes  from. 
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TABLE  5-4.  SUMMARY  OF  CONCENTRATION  DATA 


Test -Run  No.  LNG-18  (Wind  Directions  -  214°)* 


China  Lake 

Position 

Naval  Me a- 

Maximum 

Approximate 

pons  Grid 

Peak 

Time  of 

Time  of 

Time  of 

Test  Posi- 

X 

y 

z 

Concentration 

Arrival 

Peak 

Passage 

tion  No. 

(>) 

(») 

(®) 

(*) 

(sec) 

(sec) 

(sec) 

22. S 

0 

0 

48.0 

5.7 

38.3 

68.9 

32.7 

0 

0 

30.4 

9.6 

38.3 

53.6 

45.5 

0 

0 

20.6 

7.7 

26.8 

65.1 

58.2 

0 

0 

17.3 

7.7 

34.5 

S7.5 

71.0 

0 

0 

11.6 

11.5 

36.4 

76.6 

83.7 

0 

0 

9.5 

11.5 

34.4 

70.9 

96.5 

0 

0 

8.6 

11.5 

42.1 

82.3 

96.5 

44.2 

1.5 

0.0 

0 

0 

0 

105.0 

-5.1 

1.5 

5.5 

15.3 

19.2 

47.9 

109.2 

0 

0 

5.9 

11.5 

30.6 

70.9 

122.0 

0 

0 

6.9 

15.3 

47.9 

76.6 

134.7 

0 

0 

4.7 

13.4 

42.1 

76.6 

1 

47.9 

11.1 

1.5 

0.0 

9.6 

30.6 

53.6 

2 

47.9 

11.1 

1.5 

5.9 

11.5 

21.1 

53.6 

3 

75.3 

17.0 

1.5 

0.0 

11.5 

32.5 

68.9 

4 

7S.3 

17.0 

1.5 

4.0 

11.5 

38.3 

55. S 

S 

99.9 

20.4 

0.5 

3.6 

19.2 

46.0 

65.1 

6 

124.5 

25.5 

1.5 

1.7 

17.2 

30.6 

76.6 

7 

54.0 

-IS. 3 

1.5 

7.6 

13.4 

34.5 

49.8 

8 

80.3 

-53.6 

1.5 

6.4 

15.3 

38.3 

76.6 

9 

4S.S 

34.9 

1.5 

0.0 

0 

0 

0 

10 

71.0 

39.1 

l.S 

0.0 

0 

0 

0 

Test-Run 

No.  LNG-19 

(Wind 

Direction 

=  260°)* 

China  Lake 

Position 

Naval  Wea- 

Maximum 

Approximate 

pons  Grid 

Peak 

Time  of 

Time  of 

Time  of 

Test  Posi- 

X 

y 

i  Concentration 

Arrival 

Peak 

Passage 

tion  No. 

(») 

(■) 

(») 

(*) 

(sec) 

(sec) 

(sec) 

22.5 

0 

0 

48.0 

3.8 

34.5 

49.8 

32.7 

0 

0 

34.6 

7.7 

23.0 

65.1 

45.5 

0 

0 

27.2 

7.7 

26.8 

57.5 

S8.2 

0 

0 

24.3 

11.5 

34.5 

65.1 

71.0 

0 

0 

15.8 

11.5 

42.1 

88.1 

83.7 

0 

0 

15.2 

15.3 

26.8 

72.8 

96.5 

0 

0 

13.2 

19.2 

42.1 

84.3 

98.2 

-36.6 

1.5 

0.0 

0 

0 

0 

109.2 

0 

0 

11.4 

11.5 

38.3 

72.8 

122.0 

0 

0 

11.4 

15.3 

38.3 

72.8 

134.7 

0 

0 

7.8 

15.3 

38.3 

76.6 

2 

42.1 

-28.9 

0.5 

0.0 

0 

0 

0 

9 

S5.7 

-8.5 

l.S 

20.3 

7.7 

34.5 

65.1 

10 

77.8 

-23.0 

1 . 3 

0.0 

0 

0 

0 

•Coordinate  System  is  right  handed  and  the  positive  x-axis  is  orientated 
in  the  direction  of  the  wind.  Wind  Direction  is  reference  to  direction 
wind  comes  from. 
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TABLE  5-4.  (continued) 


Test-Run  No.  LNG-20  (Wind  Direction  =  256°)* 


Chin*  Lake 

Naval  Wea¬ 
pons  Grid 

Test  Posi-  x 
tion  No .  (B) 

Position 

y 

(■) 

z 

(m) 

Maxima 

Peak 

Concentration 

m 

Time  of 
Arrival 
(sec) 

Time  of 
Peak 

(sec) 

Approx iaate 
Time  of 
Passage 

(sec) 

22. S 

0 

0 

18.8 

3.8 

11.5 

70.9 

22. S 

8.5 

0 

3.9 

7.7 

19.2 

68.9 

22.5 

17.0 

0 

0.0 

0 

0 

0 

22.5 

-8.5 

0 

9.7 

5.8 

24.9 

68.9 

22.5 

-17.0 

0 

0.0 

0 

0 

0 

32.7 

0 

0 

13.2 

3.8 

23.0 

72.8 

32.7 

8.5 

0 

4.4 

7.7 

34.5 

84.3 

32.7 

-8.5 

0 

8.6 

S .  8 

42.1 

80.4 

32.7 

-17.0 

0 

0.0 

0 

0 

0 

45.5 

0 

0 

9.5 

5.8 

26.8 

84.3 

45.5 

8.5 

0 

5.6 

5.8 

11.5 

67.0 

45.  S 

-8.5 

0 

5.6 

7.7 

34.5 

72.8 

45.5 

-17.0 

0 

0.0 

0 

0 

0 

58.2 

0 

0 

7.8 

7.7 

36.4 

76.6 

58.2 

8.5 

0 

5.9 

5.8 

26.8 

72.8 

58.2 

17.0 

0 

1.1 

13.4 

23.0 

51.7 

58.2 

-8.5 

0 

5.1 

11.5 

42.1 

68.9 

58.2 

-17.0 

0 

0.5 

28.7 

28.7 

28.7 

71.0 

0 

0 

6.4 

7.7 

38.3 

72.8 

71.0 

8.5 

0 

5.4 

7.7 

38.3 

76.6 

71.0 

17.0 

0 

1.6 

13.4 

40.2 

57.5 

71.0 

25.5 

0 

0.5 

34.5 

38.3 

49.8 

71.0 

-8.5 

0 

3.6 

15.3 

30.6 

32.6 

71.0 

-17.0 

0 

0.0 

0 

0 

0 

83.7 

0 

0 

4.4 

7.7 

38.3 

76.6 

83.7 

8.5 

0 

3.9 

5.8 

26.8 

111.1 

83.7 

17.0 

0 

1.8 

11.5 

23.0 

65.1 

83.7 

25.5 

0 

0.5 

23.0 

26.8 

30.6 

83.7 

-8.5 

0 

2.9 

11.5 

21.1 

74.7 

83.7 

-17.0 

0 

0.5 

23.0 

24.9 

26.8 

96.5 

0 

0 

4.4 

7.7 

42.1 

91.9 

96.5 

8.5 

0 

3.6 

7.7 

46.0 

91.9 

96.5 

17.0 

0 

1.6 

11.5 

46.0 

65.1 

96.5 

2S.5 

0 

0.3 

95.8 

103.4 

122.6 

96.5 

-8.5 

0 

3.1 

11.5 

38.3 

80.4 

96.5 

-17.0 

0 

0.0 

0 

0 

0 

100.7 

-28.9 

1.5 

0.0 

0 

0 

0 

109.2 

0 

0 

3.4 

7.7 

42.1 

91.9 

109.2 

8.5 

0 

3.3 

7.7 

38.3 

88.1 

109.2 

17.0 

0 

1.6 

11.5 

46.0 

68.9 

109.2 

25.5 

0 

0.8 

26.8 

30.6 

34.5 

109.2 

-B.5 

0 

2.6 

7.7 

30.6 

68.9 

109.2 

-17.0 

0 

0.3 

38.3 

42.1 

42.1 

122.0 

0 

0 

2.9 

7.7 

38.3 

80.4 

134.7 

0 

0 

2.1 

7.7 

46.0 

95.8 

2  43.8 

-25.5 

0.5 

0.0 

0 

0 

0 

4  64.2 

-39.1 

1.5 

0.0 

0 

0 

0 

9  56.5 

-4.3 

1.5 

8.8 

5.8 

40.2 

99.6 

10  86.3 

-15.3 

1.5 

2.5 

3.8 

23.0 

26.8 

‘Coordinate  System  is  right  handed  and  the  positive  x-axis  is  orientated 
in  the  direction  of  the  wind.  Wind  Direction  is  reference  to  direction 
wind  comes  from. 
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TABLE  5-4. (continued) 

Test-Run  No.  LNG-21  (Wind  Direction  =  224°)* 


Chins  Lake 

Position 

Naval  Nea- 

Maximum 

Approximate 

pons  Grid 

Peak 

Time  of 

Time  of 

Time  of 

Test  Posi- 

X 

y 

z 

Concentration 

Arrival 

Peak 

Passage 

tion  No. 

(») 

(m) 

(m) 

(') 

(sec) 

(sec) 

(sec) 

122.0 

0 

0 

6.4 

15.3 

30.6 

72.8 

122.0 

8.5 

0 

5.4 

23.0 

34.5 

72.8 

122.0 

17.5 

0 

6.1 

34.5 

42.1 

57.5 

122.0 

25.5 

0 

0.8 

143.6 

153.2 

176.2 

122.0 

-8.5 

0 

8.3 

19.2 

47.9 

99.6 

122.0 

-17.0 

0 

8.8 

32.6 

55.5 

95.8 

122.0 

-25.0 

0 

8.2 

21.1 

42.1 

90.0 

122.0 

-34.0 

0 

2.1 

139.8 

153.2 

164.7 

134.7 

0 

0 

8.2 

23.0 

46.0 

88.1 

134.7 

8.S 

0 

3.6 

30.6 

51.7 

86.2 

134.7 

17.0 

0 

2.7 

38.3 

49.8 

76.6 

134.7 

25.5 

0 

0.5 

51.7 

53.6 

68.9 

134.7 

-8.5 

0 

8.6 

19.2 

57.5 

134.1 

134.7 

-17.0 

0 

8.8 

19.2 

51.7 

124.5 

134.7 

-25.5 

0 

6.4 

46.0 

42.1 

88.1 

134.7 

-34.0 

0 

2.1 

38.3 

48.0 

76.6 

1 

SO. 6 

1.7 

1. 5 

19.0 

9.6 

32.6 

72.8 

2 

50.6 

1.7 

O.S 

24.7 

11.5 

40.2 

118.7 

3 

77.8 

1.7 

1.5 

10.5 

17.2 

30.6 

88.1 

4 

77.8 

1.7 

0.5 

17.7 

13.4 

30.6 

99.6 

S 

102.4 

2.6 

O.S 

12.3 

19.2 

38.3 

91.9 

6 

127.1 

2.6 

1.5 

6.1 

15.3 

34.5 

86.2 

7 

51.4 

-23.8 

1.5 

4.3 

24.9 

28.7 

36.4 

8 

77.8 

-23.8 

1.5 

10.8 

19.2 

28.7 

68.9 

9 

SO. 6 

27.2 

1.5 

4.0 

9.6 

42.1 

67.0 

10 

76.9 

27.2 

1.5 

1.5 

0 

0 

0 

•Coordinate  System  is  right  handed  and  the  positive  x-axis  is  orientated 
in  the  direction  of  the  wind.  Wind  Direction  is  reference  to  direction 
wind  comes  from. 
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TABLE  6.  SUMMARY  OF  APPROACH  FLOW  CHARACTERISTICS 


DESCRIPTION 


zo  (m) 


1/n 


A  (m)  §  2  meters 
x 


ATMOSPHERIC  DATA* 
0.01  -  0.15 

0.143  -  0.167 
12.0  -  30.0 


MODELED  VALUES** 
0.017 

0.18 

14.5 


Az  (m)  @  2  meters 


1 


5.1 


*Counihan,  1975 

**At  wind  speed  reference  location  (see  Figure  6) 
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TABLE  7.  SUMMARY  OF  PEAK  CONCENTRATION  DATA  AT  TEST  POINT  LOCATIONS 
FOR  MODEL  AND  FIELD 


Location 
(China  Lake 


Naval  Weapons 
Grid) 

LNG- 

Field 

18 

Model 

LNG- 

Field 

19 

Model 

LNG- 

Field 

20 

Model 

LNG- 

Field 

21 

Model 

1 

>55 

5.9 

>55 

0 

1.65 

0 

1.65 

19 

2 

>S5 

5.3 

>55 

0.0 

1.0 

0.0 

1.65 

24.7 

J 

>55 

4.1 

0.75 

0 

0 

0.0 

0.7 

10.6 

4 

>5% 

4.0 

>55 

0 

0 

0.0 

1.65 

17.7 

5 

0.7 

3.6 

0 

0 

0 

0.0 

0.6 

12.3 

6 

0 
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igure  4.  Meteorological  Wind  Tunne 
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Figure  8.  Model  and  Prototype  Source  Conditions  for  LNG-21 
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Figure  10.  Velocity  Data  Reduction  Flow  Chart 
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Figure  11.  Typical  Response  of  Hot  Film  Aspirating  Probe 


Figure  12.  Hot  Film  Aspirating  Concentration  Probe 
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Figure  16-1.  Normalized  Concentration  Versus  Downwind  Distance 
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Figure 
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APPENDIX  A  -  The  Calculation  of  Model  Scale  Factors 

As  discussed  previously  in  Section  2.3  the  dominant  scaling 
criterion  for  the  simulation  of  LNG  vapor  cloud  physics  are  the 
Froude  number  and  the  volume  flux  ratio.  By  setting  these  parameters 
equal  for  model  and  prototype  one  obtains  the  following  relationships 
for  a  model  (length  scale  (L.S.)  of  1:85  and  a  model  specific  gravity 
(S.G.)  of  1.38 


In  addition  to  these  scaling  parameters  which  govern  the  flow 
physics  one  must  also  scale  the  mole  fractions  (concentrations) 
measured  in  the  model  to  those  that  would  occur  in  the  prototype. 

This  scaling  is  required  since  the  number  of  moles  being  released 
in  a  thermal  plume  are  different  than  the  number  of  moles  being 
released  in  an  isothermal  plume.  To  be  more  precise  the  relationship 
between  the  molal  flow  rate  of  source  gas  in  the  model  and  the 
propotype  is 
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n  =  (T  /T  ) 


@  b.o. 


n  =  (2.70)  n 
m  m 


By  definition  the  concentration  of  LNG  vapor  is  expressed  as: 
xp  =  nNG^nNG  +  na^ 


Substituting  model  equivalents  into  the  above  expression  yields 

(T  /T  )  n.  n, 

m  P  @b.o.  Ar  Ar 

*P  "  nTV  *  11  (T  /T  ) „ 


@b .  o . 


Ar  a  Ar  av  p  nr@b.o. 


XP  "  ».  ‘  <>  -  X.H0.37) 

This  equation  was  used  to  correct  the  modeled  measurements  to 
those  that  would  be  observed  in  the  field. 


